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Abstract. The started LHC heavy ion program makes it possible to probe new frontiers of the high temper- 
ature Quantum Chromodynamics. It is expected that the role of hard and semi-hard particle production 
processes may be significant at ultra-high energies even for bulk properties of the created matter. In 
this paper, the LHC data on multiplicity, hadron spectra, elliptic flow and femtoscopic correlations from 
PbPb collisions at center-of-mass energy 2.76 TeV per nucleon pair are analyzed in the framework of the 
HYDJETH — h model. The influence of the jet production mechanism on these observables is discussed. 
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1 Introduction 

One of the main aims of modern high energy physics is to 
study the strong interaction in extreme regimes at super- 
high densities and temperatures. Hot and dense multi- 
parton and multi-hadron systems are produced in high- 
energy nuclear collisions and can be investigated within 
the Quantum Chromodynamics (QCD) (see, e.g., recent 
reviews [1,2,3,4 ). The experimental and phenomenolog- 
ical study of multi-particle production in ultrarelativistic 
heavy ion collisions is expected to provide valuable infor- 
mation on the dynamical behavior of QCD matter in the 
form of the quark-gluon matter (QGM), as predicted by 
lattice calculations. The current lead-lead (PbPb) collision 
center-of-mass energy per nucleon pair at LHC, = 
2.76 TeV, is a factor of ~ 14 larger than that at RHIC 
and makes it possible to probe new frontiers of the super- 
high temperature and (almost) net-baryon free QCD. It is 
expected that a role of hard and semi-hard particle pro- 
duction may be significant at ultra-high energies even for 
bulk properties of the created matter. 

A number of striking LHC results from PbPb runs have 
been published by ALICE [SUSlEUHliniIinUTTlII21II3UI31IISl 
QM7], ATLAS and CMS [2~gp5P5P6| 

[27|I2^P^I^^T1[8'21[8'81[M1|M1PH] collaborations (see [37] for 
the overview of the results from the first year of heavy 
ion physics at LHC). In the paper, we analyze the LHC 
data on multiplicity, hadron spectra, elliptic flow and fem- 
toscopic momentum correlations from PbPb collisions at 
^ NN = 2.76 TeV in the framework of the HYDJET++ 
model [35] . The main goal of this study is to show a signif- 
icant influence of the (in-medium) jet production mecha- 
nism on the global event pattern at LHC. Note that the 
influence of the interactions between jets and soft hadrons 
on a number of physical observables at the LHC was stud- 
ied recently in the framework of the EPOS model [39 . 



2 HYDJET++ model 

HYDJET++ (the successor of HYDJET [40] ) is the model 
of relativistic heavy ion collisions, which includes two in- 
dependent components: the soft state (hydro-type) and 
the hard state resulting from the in-medium multi-parton 
fragmentation. The details of the used physics model and 
simulation procedure can be found in the HYDJET++ 
manual [35]. Below main features of the model are de- 
scribed only very briefly. 



2.1 Soft component 

The soft component of an event in HYDJET++ is the 
"thermal" hadronic state generated on the chemical and 
thermal freeze-out hypersurfaces obtained from the pa- 
rametrization of relativistic hydrodynamics with preset 
freeze-out conditions (the adapted event generator FAST 
MC [4"T1|42| ) . Hadron multiplicities are calculated using 
the effective thermal volume approximation and Poisson 
multiplicity distribution around its mean value, which is 
supposed to be proportional to a number of participat- 
ing nucleons for a given impact parameter of a AA colli- 
sion. To simulate the elliptic flow effect, the hydro-inspired 
parametrization is implemented for the momentum and 
spatial anisotropy of a soft hadron emission source |38l 

H3J. 

The soft hadron simulation procedure includes the fol- 
lowing steps. 

— Generation of a hadron 4-momentum in the liquid el- 
ement rest frame in accordance with the equilibrium 
distribution function. 

— Generation of a spatial position and local 4- velocity of 
the liquid element in accordance with the phase space 
and the character of fluid motion. 
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— Standard von Neumann rejection/acceptance proce- 
dure to account for a difference between the true and 
generated probabilities. 

— Boost of the hadron 4-momentum in the event center- 
of-mass frame. 

— Taking into account two- and three-body decays of res- 
onances with branching ratios taken from the SHARE 
particle decay table |44) . 

The model has a number of input parameters for the soft 
component. They are tuned from fitting to experimen- 
tal data values for various physical observables (see the 
next section) . Parameters which determine the hadron ra- 
tios are the chemical freeze-out temperature T c h and ther- 
modynamical potentials (baryon potential /ib, strangeness 
potential /is, isospin potential /xq, and potential of posi- 
tively charged pions ^ at thermal freeze-out). We have 
used T c h = 165 MeV and, for simplicity, zero thermo- 
dynamical potentials. The collective flow parameters are 
the maximal longitudinal and transverse flow rapidities, 
ymax = 4 5 and ymax = 12 q 5> respectively. The lat- 
ter value has been fixed from the data on hadron p^- 
spectra together with the thermal freeze-out temperature 
T t h — 105 MeV. The total multiplicity of the soft compo- 
nent is determined by space-time parameters on the ther- 
mal freeze-out for central collisions, maximal transverse 
radius Ri(b = 0), proper time Tf(6 = 0) and the dura- 
tion of hadron emission Arf(b = 0). These parameters 
have been fixed from fitting of three-dimensional correla- 
tion functions measured for tt + tt + pairs and extracting the 
correlation radii: R { (b = 0) = f3.45 fm, r f (6 = 0) = 12.2 
fm/c, Ar f (6 = 0) = 2 fm/c. 

2.2 Hard component 

The model used for the hard component in HYDJET++ is 
the same as in the HYDJET event generator. It is based 
on the PYQUEN partonic energy loss model [40]. The 
approach describing the multiple scattering of hard par- 
tons is based on accumulated energy loss via gluon ra- 
diation which is associated with each parton scattering 
in expanding quark-gluon fluid. It also includes the in- 
terference effect in gluon emission with a finite forma- 
tion time using the modified radiation spectrum dE/dx 
as a function of the decreasing temperature T. The model 
takes into account radiative and collisional energy loss 
of hard partons in longitudinally expanding quark-gluon 
fluid, as well as the realistic nuclear geometry. Radiative 
energy loss is treated in the framework of the BDMS 
model 45,46 with a simple generalization to a massive 
quark case using the "dead-cone" approximation [UJ . The 
collisional energy loss due to elastic scatterings is cal- 
culated in the high-momentum transfer limit [151149 [I5U] . 
The strength of the energy loss in PYQUEN is deter- 
mined mainly by the initial maximal temperature Xj!j nax 
of hot matter in central PbPb collisions (i.e. an initial 
temperature in the center of the nuclear overlap at mid- 
rapidity). The energy loss depends also on the proper time 
tq of QGM formation and the number N{ of active flavors 



in the medium. The parameter values T™ ayi = 1 GeV, 
To = 0.1 fm/c and Nf = (gluon-dominated plasma) are 
used at \/*nn = 2.76 TeV. Another important ingredient 
of PYQUEN is the angular spectrum of in-medium radia- 
tion. Since the ATLAS [TSJ and CMS [H] data on the dijet 
asymmetry support the supposition of the "out-of-cone" 
medium-induced partonic energy loss, the "wide-angular" 
parametrization of medium-induced gluon radiation |51j 
is used in the presented investigation. 

The simulation of single hard nucleon-nucleon sub- 
collisions by PYQUEN is constructed as a modification 
of the jet event obtained with the generator of hadron- 
hadron interactions PYTHIA_6.4 [52]. We used PYTHIA 
tune Pro-Q20 for the present simulation. It reproduces 
the CMS data on hadron momentum spectra in pp colli- 
sions at LHC energies with the 10-15% accuracy in the 
full measured px-range JS2| ■ 

The event-by-event simulation procedure in PYQUEN 
includes the following steps. 

— Generation of initial parton spectra with PYTHIA and 
production vertices at a given impact parameter. 

— Rescattering-by-rescattering simulation of the parton 
path in a dense zone and its radiative and collisional 
energy loss. 

— Final hadronization according to Lund string model 
for hard partons and in-medium emitted gluons. 

Then the number of PYQUEN jets is generated accord- 
ing to the binomial distribution. The mean number of jets 
produced in an AA event is calculated as a product of the 
number of binary NN sub-collisions at a given impact pa- 
rameter per the integral cross section of the hard process 
in NN collisions with the minimum transverse momen- 
tum transfer p™ n . The latter is an input parameter of the 
model. In the framework of HYDJET, partons produced in 
(scmi)hard processes with the momentum transfer lower 
thanp™ 1 ", are considered as being "thermalized" . So, their 
hadronization products are included "automatically" in 
the soft component of the event. In order to take into 
account the effect of nuclear shadowing on parton distri- 
bution functions, we use the impact parameter dependent 
parametrization obtained in the framework of Glauber- 
Gribov theory |54j . The nuclear shadowing correction fac- 
tor S(rx, ^2, Xi, X2, Q 2 ) depends on kinematic variables of 
incoming hard partons (the momentum fractions x\ % of 
initial partons from incoming nuclei and the momentum 
scale Q 2 of a hard scattering), and on the transverse coor- 
dinates rx t 2 of partons in their parent nuclei. This initial 
state effect reduces the number of partons in the incoming 
hadronic wave-functions of both nuclei and, thus, reduces 
the yield of high-pT hadrons coming from jets (in addition 
to the medium-induced partonic energy loss). 

To study the femtoscopic (HBT) momentum correla- 
tions, it is necessary to specify a space-time structure of a 
hadron emission source. The information about final par- 
ticle 4-momenta and 4-coordinates of emission points is 
automatically recorded using the soft hadron simulation 
procedure (see previous subsection). For the hard compo- 
nent, the existence of a "source" with a given size is not 
obvious and is a matter of discussion (see, for instance, 
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[55] ). The point of jet hadronization and the point of ini- 
tial parton-parton hard scattering do not coincide. In ac- 
cordance with [55] . we assume that hadronization occurs 
at different distances from an initial hard scattering, de- 
pending on an energy of a jet. 

Our treatment of the coordinate information for low 
momentum jet particles with px < 1 GeV/c is similar to 
one for directly produced hadrons from soft component. 
Such particles are emitted from the fireball of radius R{ 
at mean proper time Tf with the emission duration Art . 

Four-coordinates of high momentum jet particles with 
px > 1 GeV/c are shifted from those of a jet vertex and 
coded in the following way [55] : 



cos < 
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where the subscripts i and j refer to particles and jets 
respectively. Here xq = r cos ip and yo — r sin tp are trans- 
verse coordinates of a jet production vertex, r is a distance 
from the nuclear collision axis to the vertex with the az- 
imuthal angle ip. The longitudinal position and the initial 
time of a hard process are estimated as zq ~ smhyj/p^j 
and to — coshy j /pTj, respectively (<pj and yj are the 
azimuthal angle and the pseudorapidity of a jet, respec- 
tively). For each hadron i, the localization Lu of hadroniza- 
tion along the direction of a jet is randomized from Gaus- 
sian distributions with mean lj and variance a — [55] . 
The hadronization points L t u and L t 2i (with respect to 
the jet direction) are randomized in the transverse di- 
rection from a Gaussian distribution with variance cr t = 
0.5 fm and zero mean value; U is assumed to be propor- 
tional to jet momentum |p -|. Thus, U = f\pj\ where / is 
a multiplicative factor that represents our lack of theoret- 
ical insight into the process of hadronization (the value 
/ = 5 being used in our simulations). For simplicity, the 
jet axis is determined by the momentum direction of a 
leading particle in the event. Then four-coordinates of jet 
particles with px > 1 GeV/c are shifted from the jet ver- 
tex position along or opposite this direction depending on 
a sign of the scalar production (pipj). 



3 Numerical results and discussion 

It was demonstrated in [35] that the HYDJET++ model 
can describe bulk properties of the hadronic state created 
in AuAu collisions at RHIC at y/s NN — 200 GeV (hadron 
spectra and ratios, radial and elliptic flow, femtoscopic 
momentum correlations), as well as the high-px hadron 
spectra. At present we apply HYDJET++ with tuned in- 
put parameters to reproduce the LHC data from PbPb 
collisions, and to estimate an influence of the hard pro- 
duction mechanism on physics observables. 



3.1 Centrality and pseudorapidity dependence of 
multiplicity 

Figure [1] shows the charged multiplicity density at mid- 
rapidity (normalized on the mean number of participat- 
ing nucleons (iV par t)) as a function of event centrality 
(left) and pseudorapidity distribution in two centrality 
bins (right) in PbPb collisions at yfs NN = 2.76 TeV. The 
results of HYDJET++ simulation are compared with the 
ALICE [7J and CMS [27] data. Since multiplicities from 
soft and hard components depend on the centrality in a 
different way, the relative contribution of soft and hard 
components to the total event multiplicity can be fixed 
through the centrality dependence of dN/drj. The latter 
can be reproduced at ~ 25% contribution of the hard com- 
ponent (versus ~ 15% at RHIC [35]) in most central col- 
lisions, and it decreases for more peripheral events. This 
contribution of the hard component corresponds to the 
minimal transverse momentum transfer of parton-parton 
scatterings p™ 111 = 8.2 GeV/c. The approximately flat 77 
dependence of the multiplicity observed by CMS (the vari- 
ation being less than 10% in the pseudorapidity range 
\r/\ < 2.5) is also well reproduced by HYDJET++ sim- 
ulation for different event centralities at reasonable values 
of maximal longitudinal flow rapidity, y^ lax = 3.5 + 4.5. 

The ratios of mid-rapidity yields of final charged had- 
rons in the model (excluding weak decays of strange had- 
rons) are if ± /7r ± = 0.153 and p ± /n ± = 0.065 at 0-5% of 
PbPb centrality. The calculated kaon-to-pion ratio repro- 
duces the experimental number reported by ALICE [56 
(0.155 ±0.012), while the proton-to-pion ratio is overes- 
timated in the model by ~ 40% as compared with the 
ALICE data (0.046 ± 0.004). Since pure thermal models 
meet with difficulties to reproduce these hadron ratios, in 
our opinion, it may be interpreted as a significant influ- 
ence of (mini-)jct production on a particle composition 
in heavy ion collisions (hadron ratios for soft and hard 
components are quite different). 



3.2 Transverse momentum spectrum and nuclear 
modification factor 

Figure [2] (left) shows the transverse momentum spectrum 
of charged particles (\rj\ < 0.8) in 5% of most central 
PbPb collisions at v^nn = 2 - 76 TeV - The results of HYD- 
JET++ simulation are compared with the ALICE data [8j. 
The slope of px-spectra allows one to fix the thermal 
freeze-out temperature T th = 105 MeV and the maxi- 
mal transverse flow rapidity Y™ ax = 1.265. One can see 
that HYDJET++ reproduces the measured px-spectrum 
for central PbPb events at least up to eight orders of cross 
section magnitude. Multiplicities of the hard and soft com- 
ponents become comparable at px ~ 4.5 GeV/c. 

High transverse momentum region of hadron spectra 
(px > 5-^10 GeV/c) is sensitive to parton production and 
jet quenching. In HYDJETH — (-, they are determined by 
the PYQUEN parameters for initial conditions of QGM 
formation [51) . The nuclear modification factor Raa is 
defined as a ratio of particle yields in AA and pp collisions 
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normalized on the mean number of binary nucleon-nucleon 
sub-collisions (JV co u): 



NX 



2.76 



Raa(pt) 



d 2 N AA /dr]dpT 
{N col i) d 2 NPP/dr)d Pl 



(2) 



In the absence of nuclear effects (in initial or final states) 
at high pt, it should be Raa = 1< In figure [2] (right), 
the calculated nuclear modification factor Raa {pt ) for 
charged hadrons (\r]\ < 1) is compared with the CMS 
data [29 in 5% of most central PbPb collisions. The rea- 
sonable agreement with the data is achieved up to pt ~ 
100 GeV/c. However the simulated px-dependence of Raa 
looks rather weaker than the tendency of the data. It could 
indicate that the real energy dependence of partonic en- 
ergy loss is somewhat weaker than it is expected from the 
model simulations. It is worth to remind here that the 
energy dependence of BDMS radiative loss [dTflHS] (the 
dominant mechanism of medium-induced energy loss in 
PYQUEN) is between AE oc \nE and AE cx \ r E (de- 
pending on the energy of the radiated gluon and on the 
properties of the medium). 



3.3 Femtoscopic momentum correlations 

Effects of quantum statistics and final state interactions 
make the momentum (HBT) correlation functions of two 
or more particles in their c.m.s. sensitive to space-time 
characteristics of production process on the level of fem- 
tometers at small relative momenta. Since HYDJET++ 
specifies the space-time structure of a hadron emission 
source (for both soft and hard components), the momen- 
tum correlation function can be introduced by the stan- 
dard weighting procedure [57,58 . Having the information 
about final particle 4-momenta pi and 4-coordinates x% of 
emission points, we can calculate the correlation function 
with the help of the weight procedure, assigning a weight 
to a given particle combination and accounting for effects 
of quantum statistics: 



1 + cos(g • Ax), 



(3) 



where q = p\ — p2 and Ax = x\—%i. Then the correlation 
function is defined as a ratio of a weighted histogram of 
pair kinematic variables to an unweighted one. The corre- 
sponding correlation widths are parametrized in terms of 
the Gaussian correlation radii Ri as follows: 



CF( Pl ,p 2 ) 



1 + Xexp(-Rl ut ql ut 



r>2 2 
-^-sidcysidc 

-^long^long 



(4) 



where q = (q ou t, Qsidet 9iong) is the relative three-momen- 
tum vector of two identical particles and A is the corre- 
lation strength. The "out" axis is pointing along the pair 
transverse momentum, the "side" axis is perpendicular to 
the "out" direction in the transverse plane, and the "long" 
axis is along the beam. 

Figure [3] shows the fitted correlation radii as a function 
of mean transverse momentum fcx of pion pair (|?/| < 0.8) 



in 5% of most central PbPb collisions at y/s 
TeV. The results of the HYDJET++ simulation are com- 
pared with the ALICE data [5] ■ Taking into account both 
soft and hard components allows us to reproduce the mea- 
sured fcx-dependencies of the correlation radii. 

Let us comment how our relatively simple model with- 
out detailed fluid dynamic evolution and hadronic cas- 
cade is able to fit the experimental data on the femto- 
scopic observables. In fact, there are 3D-hydrodynamics 
models which reproduce HBT radius parameters in heavy 
ion collisions at RHIC and LHC energies, e.g. [ST)II60| . 
A successful simultaneous description of hadronic yields 
and spectra, elliptic flow and correlation radii at RHIC 
and LHC is achieved also in hydro- kinetic model [BTj . 
These approaches include the hydrodynamical evolution 
and hadronic cascade calculations which are very com- 
puter time consuming. HYDJETH — h is the multi-purpose 
event generator for a fast simulation (needed e.g. in Monte- 
Carlo studies of the experimental setup), and so it does not 
contain detailed hydrodynamical evolution and hadronic 
cascade. The treatment of the soft component in HYD- 
JET++ is based on the fast simulation of the freeze-out 
hypersurface realized in FAST MC model 41,42 describ- 
ing well the femtoscopic radii at RHIC energies. The main 
feature of this approach is that the hadron system expands 
hydrodynamically with the frozen chemical composition, 
cools down, and finally decays at some thermal freeze- 
out hypersurface. The chemical potentials are fixed from 
the experimental particle number ratios at the chemical 
freeze-out. In spite of lack of the hydrodynamical evo- 
lution in HYDJET++, the assumption of the conserva- 
tion of the particle number ratios from the chemical to 
thermal freeze-out allows us to calculate the chemical po- 
tentials at the thermal freeze-out. The absolute values of 
particle densities are determined by the free parameter 
of the model, effective pion chemical potential ^ at the 
thermal freeze-out. The best tune of the temperatures at 
chemical and thermal freeze-out hypersurfaces and vol- 
ume parameters allow us to reproduce simultaneously the 
data on the total multiplicity, hadronic spectra and cor- 
relation radii. At the same time, low momentum particles 
(with pt < 1 GeV/c) coming from the hard component 
are emitted from the same thermal freeze-out hypersurface 
as soft component hadrons (see subsection 2.2). 



3.4 Elliptic flow 

The elliptic flow coefficient v 2 is defined as the second- 
order Fourier coefficient in the hadron distribution over 
the azimuthal angle if relative to the reaction plane an- 
gle V'R) so that V2 = (cos2(y;— V'r))- It is an important 
characteristic of physics dynamics at early stages of non- 
central heavy ion collisions. According to the typical hy- 
drodynamic scenario, values of v 2 (pt) at low-px are de- 
termined mainly by internal pressure gradients of an ex- 
panding fireball during the initial high density phase of 
the reaction, while the elliptic flow at high-px is generated 
due to partonic energy loss in an azimuthally asymmetric 
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volume of QGM. Figures H] and [5] show elliptic flow coeffi- 
cient V2 as a function of the hadron transverse momentum 
Pt and pseudorapidity 77, respectively. The results of the 
HYDJET++ simulation are compared with the CMS data 
obtained from the 4-particle cumulant method [33J . HYD- 
JET++ reproduces ^-dependence and pr dependence of 
V2 at mid-rapidity (up to px ~ 5 GeV/c and 40% cen- 
trality). As it is expected, taking into account the hard 
component changes the hydrodynamical growth of V2 with 
Pt, so V2 decreases at pt > 3^4 GeV/c. The flat in- 
dependence of i>2 manifests the strong longitudinal flow 
effect taken into account in HYDJET++. 



4 Conclusions 

The LHC data on multiplicity, hadron spectra, elliptic flow 
and femtoscopic momentum correlations in PbPb colli- 
sions at yi NN = 2.76 TeV are analyzed in the frame- 
work of the HYDJET++ model. The significant influence 
of the jet production mechanism on these observables is 
shown. We have estimated that the contribution of the 
hard component to the total multiplicity is about 25% at 
mid-rapidity in most central collisions, and it decreases 
for more peripheral events. 

Taking into account both hard and soft components 
and tuning input parameters allow HYDJET-I — h to re- 
produce the LHC data on centrality and pseudorapidity 
dependence of charged particle multiplicity, px-spectra, 
nuclear modification factor (up to pt ~ 100 GeV/c) and 
n^ir^ correlation radii in central PbPb collisions, and pt- 
and ^-dependencies of the elliptic flow coefficient V2 (up 
to pt ~ 5 GeV/c and 40% centrality). 

Discussions with I. Arsene, L. Bravina, D. d'Enterria, A. Gribu- 
shin, V. Korotkikh, L. Sarycheva, Yu. Sinyukov, K. Tiwonyuk 
and E. Zabrodin are gratefully acknowledged. We thank col- 
leagues from CMS and ALICE collaborations for fruitful coop- 
eration. This work was supported by Russian Foundation for 
Basic Research (grants 10-02-93118 and 12-02-91505), Grant 
of President of Russian Federation for scientific Schools sup- 
porting (3920.2012.2) and Dynasty Foundation. 
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Fig. 1. Charged multiplicity density at mid-rapidity (normalized on the mean number of participating nucleons (iVpart}) as a 
function of PbPb event centrality (left) and pseudorapidity distribution in two centrality bins (right) at \/s NN = 2.76 TeV. The 
open and closed points are ALICE [7] and CMS 27! data respectively, curves (left) and histograms (right) are the simulated 
HYDJET++ events (solid - total result, dashed - hard component, dotted - soft component). 




10 12 14 16 



p , GeV/c 




10* 

p , GeV/c 



Fig. 2. Charge particle transverse momentum spectrum (left) and nuclear modification factor Raa(pt) (right) in 5% of most 
central PbPb collisions at v^nn = 2.76 TeV. The points are ALICE [8] (left) and CMS |29| (right) data, histograms are the 
simulated HYDJETH — h events (solid - total result, dashed - hard component, dotted - soft component). 
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Fig. 3. tt^tt^ correlation radii as functions of pion pair transverse momentum ki in 5% of most central PbPb collisions at 
V^nn = 2.76 TeV. The points are the ALICE data [5], curves are the simulated HYDJET++ events. 
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Fig. 4. Elliptic flow coefficient V2(pt) for charged hadrons at mid-rapidity (\r)\ < 0.8) for different centralities of PbPb collisions 
at %/®nn = 2-76 TeV. The points are the CMS data [33], histograms are the simulated HYDJETH — h events (solid - total result, 
dashed - hard component, dotted - soft component). 
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Fig. 5. Elliptic flow coefficient V2(rj) (0.3 < pr < 3 GeV/c) for charged hadrons at different centralities of PbPb collisions at 
y/s NN = 2.76 TeV. The points are the CMS data [33], histograms are the simulated HYDJET++ events. 



